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Abstract. Red supergiants with their enormous brightness at J-band 
are ideal probes of cosmic chemical composition. It is therefore cru- 
cial to have realistic models of radiative transfer in their atmospheres, 
which will permit determination of abundances accurate to 0.15 dex, 
the precision attainable with future telescope facilities in galaxies as 
distant as tens of Mpc. Here, we study the effects of non-local ther- 
modynamic equilibrium (NLTE) on the formation of iron, titanium, 
and silicon lines, which dominate J-band spectra of red supergiants. 
It is shown that the NLTE radiative transfer models enable accurate 
derivation of metallicity and effective temperature in the J-band. We 
also discuss consequences for RSG spectrum synthesis in different spec- 
tral windows, including the heavily TiO-blanketed optical region, and 
atmospheric structure. We then touch upon challenges of NLTE inte- 
gration with new generation of 3D hydrodynamical RSG models and 
present the first calculations of NLTE spectra with the mean 3D model 
of Betelgeuse. 

Introduction 

m ^ . Chemical composition is one of the key observable characteristics of star forming 
' galaxies in the nearby and in the high redshift universe. So far, most of our inter- 
s' , mation about their metal content has been obtained from the analysis of strong 
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emission lines from H II regions. However, measurements of galaxy metallicities 
are then uncertain by a large factor because of the systematic uncertainties inher- 
ent to this 'strong-line' method. Furthermore, the method yields basically only the 
oxygen abundance, which is then taken as a placeholder for the overall metallicity. 
In this case, there is no information on abundance ratios, which can be a powerful 
diagnostic of the chemical enrichment history. 

An alternative approach avoiding these weaknesses is the spectroscopic anal- 
ysis of supergiant stars. Here, much progress has been made through the optical 
spectroscopy of blue supergiants in the Milky Way and few other galaxies of the 
Local group. For extragalactic astrophysics, however, red supergiants (RSGs) are 
more promising candidates. The spectral energy distribution of RSG's peaks in 
the infrared, where interstellar extinction is reduced. Particularly attractive for 
quantitative spectroscopy is the J-band, which contains many isolated atomic lines. 
Spatial resolution in the infrared is also higher than in the optical thanks to more 
efficient corrections by adaptive optics systems. Extragalactic RSGs are therefore 
ideal targets for spectroscopy with future telescope facilities, such as the European 
Extremely Large Telescope. The abundances of various chemical elements could 
be then directly measured out to distances of 70 Mpc, far beyond our Local Group 
of galaxies. 

The analysis of RSG spectra, however, is a challenging task. The major com- 
plexities arise due to very low gravities, stipulating departures from Local Thermo- 
dynamic Equilibrium (LTE) in their photospheres, and violent surface convection. 

We present a detailed investigation of non-LTE (NLTE) line formation in RSG 
atmospheres using new atomic models of Fe, Si, and Ti computed with state- 
of-the-art atomic data. The radiative transfer calculations are performed with 
ID hydrostatic and mean 3D hydrodynamic model atmospheres. The effects on 
spectrum synthesis and element abundance determinations from the atomic lines 
important in the near-IR spectra of RSG's are quantified. 

2 Line formation and atomic models 

The atomic models of Ti, Fe, and Si have been described elsewhere (Bergemann 
et al. 2012, 2 013|) . In short, they were constructed from the best available atomic 
data, including the experimental data from NIST, theoretically-predicted energy 
levels and radiative transitions from the Kurucz's database, and photoionization 
cross-sections from TOPbase. The model atoms were tested by performing test 
calculations for several benchmark stars with well-known physical parameters, 
including the Sun and Arcturus. Fig. Q] shows the Si I atom. 

As a basis for NLTE statistical equilibrium calculations, we use two types of 
model atmospheres. The MARCS model atmospheres (Gustafsson et al. 2008) are 
ID hydrostatic models, widely used in spectroscopy of cool stars. We also attempt 
the non-LTE line formation with the spatially-averaged models constructed from 
the C0 5 BOLD star-in-a-box simulations of stellar surface convection (see below). 
The model tested in this work is that from Chiavassa et al. (|20 1 1[) . st35gm03nl3, 
with luminosity 89000Lsun and T e g = 3430 K, and the <3D> model is computed 
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Fig. 1. The Grotrian diagram of the Si I model atom. 

by a spherical averaging of the original simulation (Fig. The details of the 
averaging procedure will be provided in a forthcoming publication. This is the 
first time that such models are coupled to NLTE radiative transfer with the goal 
to understand the impact on the RSG spectroscopic analysis. The DETAIL code 
is used for the calculation of NLTE occupation numbers. 

Line profiles are computed with the code SIU using the level departure coef- 
ficients from DETAIL. SIU and DETAIL share the same physics and line lists. 
For all details such as model structure and geometry, background opacities, solar 
abundance mixture, etc. we refer the reader to Bergemann et al. 1)20121 f2013|) . 



3 Abundances from non-LTE and LTE models 

In what follows, we first discuss the results obtained with the MARCS ID hydro- 
static model atmospheres because a large grid of such models for RSG parameters 
is available and systematic vartiations of NLTE effects with T e g, log g, and [Fe/H] 
can be well-understood. Given the very complex structure of the atoms, we focus 
only on the energy levels and transitions related to the near-IR J-band analysis of 
RSG stars (Si I: multiplet 4s 3 P° - 4p 3 D, Ti I: b 3 F - z 3 D°) 

The inclusion of NLTE effects changes the titanium abundances dramatically 
compared to LTE and with different abundance correction factors A for different 
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log x (500 nm) 

Fig. 2. The ID MARCS and <3D> model atmospheres used for NLTE radiative transfer. 
Stellar parameters (T e s, \ogg, and [Fe/H]) are indicated in the plot. 



temperatures and metallicities. The NLTE abundance correction^) vary smoothly 
between —0.4 dex and +0.2 dex as a function of effective temperature, which is 
caused by the interplay of continuum processes and line transition cascades to and 
from the relevant atomic energy levels. The NLTE effects on the J-band Fe I lines 
are much smaller. The abundance corrections become non- negligible, A > 0.1 dex, 
only for T c g > 4000 K, when radiative pumping to the higher levels takes over. 

In what concerns the J-band Si I lines, the situation resembles the classical case 
of the resonance line formation: photo-absorption in the line transition and com- 
peting de-population of the upper levels by spontaneous emission and collisional 
de-excitation. The NLTE abundance corrections are substantial and are always 
negative, with values between —0.4 to —0.1 dex. In contrast to Ti I and Fe I, the 
NLTE effects in the Si I J-band lines are weaker at higher effective temperature. 



lr The NLTE correction is a logarithmic correction, which has to be applied to an LTE abun- 
dance determination of a specific line to obtain the correct value corresponding to the use of 
NLTE line formation. We calculate these corrections at each point of our model grid for each 
line by matching the NLTE equivalent width through varying the Si abundance in the LTE 
calculations. 
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4 Non-LTE and mean 3D RHD model atmospheres 

At present, it is not possible to systematically evaluate the effect of the ID hy- 
drostatic equilibrium approximation on the determination of abundances from the 
observed spectra of red supergiants. For such calculations, full grids of 3D RHD 
models in the relevant space of stellar parameters are necessary, which are not yet 
available. 

The only realistic non-gray 3D radiative hydrodynamics simulation of convec- 
tion for a single RSG model was performed by Chiavassa et al. (2011). They 
demonstrated that there is no unique ID static model which would provide sim- 
ilar observable characteristics, i.e., SED shape and colors, as the analogous 3D 
RHD model. To match the optical spectrum dominated by TiO opacity, a ~ 100 
K cooler ID model is needed, whereas the IR spectral window can be best matched 
using a ~ 200 K hotter ID model. The NLTE effects, however, are sensitive to the 
radiation field across the whole spectrum, from UV to the IR. Thus, it is not pos- 
sible to predict the changes in the NLTE abundance corrections simply by using 
a cooler or hotter ID hydrostatic model. 

We have performed full NLTE spectrum synthesis with the <3D> model atmo- 
sphere of a cool RSG star (see Sec. 2). The model has T c g = 3430 K, log g = —0.34, 
and [Fe/H] = 0. For comparison, we used a MARCS ID static model computed for 
the same physical parameters albeit a lower effective surface gravity, \ogg = —0.5, 
comparable to what is necessary to account for the effects of turbulent pressure 
in the static models (Chiavassa et al. 1201 1|) . As expected, the differences between 
the <3D> and ID static results depend critically on the type of NLTE effects 
participating in the distribution of atomic level populations. This is illustrated in 
Fig. [3] and Fig. [4] for the lines of Si I and Ti I, respectively. 

Si I is a majority ion because of its very high first ionization potential, E ~ 8.15 
eV. The NLTE effects in the high-excitation (4.9 eV) Si I line are entirely due to 
the radiation field in the line itself, and variation of the background radiation field 
emerging from the deep atmosphere has no effect. The line formation takes place 
in the optical depth range — 2 < logrsoonm < 0, where the thermal and pressure 
structures of the ID and <3D> models are identical. As a consequence, the Si I 
line profiles computed with the MARCS and the <3D> are very similar and both 
become notably darker in the cores compared to LTE. This implies negative NLTE 
abundance corrections in either case. 

In contrast, for the near-IR low-excitation (1.4 eV) Ti I lines the effect is 
different. Ti I is a minority ionization stage (Nxii/iVriii ~ 10 -2 ) at all optical 
depths. Thus, photoionization and pumping by the radiation field originating 
from the atmospheric depths, where the <3D> models is ~ 200 K warmer than 
the ID model, is amplified. This causes the NLTE line weakening so that an LTE 
analysis with ID or <3D> models will strongly overestimate the abundance. 

We also performed LTE radiative transfer calculations with the <3D> models 
in the optical region centered on the TiO molecular bands. The lines form at in- 
termediate optical depths and, given the similarity of the mean model structures, 
no significant differences were found between both models in LTE. This result is 
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Fig. 3. NLTE line formation with the ID static and <3D> model atmosphere for the Si 
I line at 4.8 eV. 



not unexpected. Molecules are very sensitive to temperature and pressure inho- 
mogeneities (Uitenbroek & Criscuoli I201ip , thus the correct physical picture of 
TiO line formation demands NLTE calculations with full 3D RHD simulations. 

5 Conclusions 

We investigated the effects of non-local thermodynamic equilibrium on the forma- 
tion of strong Fe I, Ti I, and Si I lines, which are very prominent in the J-band 
spectra of red supergiants (RSG). Since our long-term goal is a comprehensive 
near-IR analysis of RSG's, only the model atmospheres with T e s < 4400 K and 
logs < 1 are considered. 

LTE models are fully inadequate for spectral analysis of red supergiants, both 
with ID hydrostatic and 3D hydrodynamic models. The NLTE radiative transfer 
models lead to significantly different line profiles compared to LTE. These differ- 
ences are most pronounced if the averaged 3D RHD models are used, which is also 
known from the <3D> NLTE analysis of warmer FGK stars (Bergemann et al. 
I2012|) . The synthetic NLTE line profiles and abundances inferred using these are 
similar in ID and <3D> for the high-excitation Si I lines, but they are radically 
different for the low-excitation Ti I lines. It is thus recommended, in order not to 
introduce systematic effects in determinations of RSG chemical composition, to 
avoid low-excitation lines of neutral metals in the near-IR. 

LTE spectrum synthesis in the optical spectral window, centered on the heavily 
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Fig. 4. NLTE line formation with the ID static and <3D> model atmosphere for the 
Ti I line at 1.4 eV. 



TiO-blanketed region, did not reveal any significant differences between using the 
ID static and the <3D> model. NLTE calculations with full 3D RHD simulations 
of RSG atmospheres are necessary in this case. 
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